Coulomb Gap in Graphene Nanoribbons 
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We investigate the density and temperature-dependent conductance of graphene nanoribbons with 
varying aspect ratio. Transport is dominated by a chain of quantum dots forming spontaneously due 
to disorder. Depending on ribbon length, electron density, and temperature, single or multiple quan- 
tum dots dominate the conductance. Between conductance resonances cotunneling transport at the 
lowest temperatures turns into activated transport at higher temperatures. The density-dependent 
activation energy resembles the Coulomb gap in a quantitative manner. Individual resonances show 
signatures of multi-level transport in some regimes, and stochastic Coulomb blockade in others. 



Monolayer graphene shows impressive material stabil- 
ity, even if shaped into nanostructures of about 10 nm in 
size [IHS]. Its electronic properties are tunable by gate 
electrodes [H [5] like conventional semiconductors, while 
its conductivity competes with that of metals. Graphene 
nanoribbons have the potential to be used in nanoelec- 
tronics [2^ , and graphene nanoconstrictions are the basic 
building blocks for quantum devices [3]. 

The transport properties of graphene ribbons and con- 
strictions on a Si02 substrate have been one of the puz- 
zles for the understanding of graphene nanostructures. 
Early predictions of an energy gap in ribbons [&- 8 have 
triggered intense experimental |9UT9] and theoretical re- 
search [SI [2QH3T] . It has become evident experimentally 
that localized states due to edge and bulk disorder sup- 
press the conduction and lead to a transport gap [11 - 
fT3] rather than a true band gap. In addition, exper- 
iments indicate the formation of an interaction driven 
Coulomb gap pTHl5] . A wealth of theoretical ideas rang- 
ing from Anderson localization [23H27] to Coulomb block- 
ade [21] 1^ try to explain the phenomenology. 

We show in this Letter that electronic transport in nar- 
row nanoribbons is dominated by a chain of one or mul- 
tiple quantum dots forming due to disorder. Not only 
is the conductance activated between conductance reso- 
nances, but the activation energy at each density corre- 
sponds to the Coulomb gap. At the lowest temperatures, 
cotunneling is present. Our experiments indicate that 
transport through graphene nanoribbons can be under- 
stood based on the mesoscopic details of the sample in 
a single-particle picture including Coulomb blockade. In 
contrast to recent suggestions [151 Hi] i there is no indi- 
cation that additional energy scales or mechanisms are 
necessary to describe the observed behavior. 

Graphene nanoribbons with widths below 120 nm 
and lengths of 100 nm and 200 nm were fabricated 
as described in Ref. [13] on a Si02 layer covering the 
highly doped Si-substrate which serves as a global back- 
gate. Five different devices [length (nm) x width (nm): 
200x75, 100x45, 100x80, 100x100, 100x120] were char- 
acterized in detail within this study, all showing the same 
qualitative behavior. 

The measurements were carried out in the variable 
temperature insert of a ^He cryostat with a base temper- 



FIG. 1: (a) Scanning force micrograph of the nanoribbon in- 
vestigated here (L = 200 nm, W = 7^ nm). (b) G'(Vbg) in 
a large density range showing the transport gap around Vbg 
= -2 V. The measurement was taken at T = 1.25 K with a 
source-drain bias of Vbias = 500 /iV. (c) Finite-bias measure- 
ment inside the transport gap (same temperature as (b)). 



ature of 1.25 K. The conductance was measured using 
standard lock-in techniques at 13 Hz. 

In the following we limit the detailed presentation 
of the results to the representative device displayed in 
the scanning force micrograph in Fig. [lla) with ribbon 
length L = 200 nm and width W = 75 nm. Changing the 
back-gate voltage from hole transport to electron trans- 
port allows us to locate the charge neutrality point to be 
around -2 V in back-gate voltage [see Fig. [ijb)]. Like in 
earlier studies [9Ul9] , a region of suppressed conductance 
is present around this gate voltage [shaded region in Fig. 
[lib)]. This regime is commonly referred to as the trans- 
port gap and gives an estimate for the amplitude of the 
potential inhomogeneity in the ribbon [lTHT3l. 

The behavior of the conductance in Fig. [lib) for 
Fbg > -2 V is qualitatively similar to earlier observa- 
tions in narrow disordered channels in Si-inversion lay- 
ers [32 , where the large conductance fluctuations at low 
charge carrier densities were attributed to structure in 
the density of states leading to hopping transport be- 
tween strongly localized states. They are smeared out 
as either the temperature or the charge carrier density 



is increased. Inside the transport gap the smah value of 
the conductance G <^ e'^/h indicates that the system is 
strongly localized [33 . In the investigated device the size 
of the gap is AI^g ^ 3.5 V in good agreement with the 
statistics from other measurements and the scaling law 
introduced by Han et al. [9 , which relates the width of 
the nanoribbon with the energy of the transport gap. The 
studies in Refs. 1161 and ITS have shown that the transport 
gap is largely independent of the ribbon length. 
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FIG. 2: (a) T-dependence of G inside the transport gap at 
Vbias = 100 /jM. Different curves are taken at T = 1.25 K to 
45 K (black to red lines). Inset: Coulomb resonances (grey 
line) reconstructed by a convolution of three Lorentzians with 
the derivative of the Fermi distribution (red dotted line) . (b) 
Zoom into two exemplary peaks of (a). The left peak is broad- 
ened and grows with increasing T and the right peak exhibits 
an overall decrease of G with temperature, (c) C as a func- 
tion of 1/T at three positions in Vbg indicated by arrows in 
(b). Solid lines are fits to the data according to Eq. (IT]). 



For further characterization of the device we measured 
the conductance inside the transport gap at finite-biases 
applied between source and drain. The recorded dia- 
monds of suppressed current are shown in Fig. file). Di- 
amonds of different sizes can be identified which some- 
times overlap. However, in some regimes (e.g. around 
^BG = -1-5 V) resonances at zero source-drain bias are 
observable which separate adjacent diamonds from each 
other indicating single quantum dot behavior rather than 
transport through multiple dots. This phenomenology is 
usually referred to as stochastic Coulomb blockade [34 . 

Coulomb interactions play an important role in 
graphene and lead to the formation of a Coulomb gap 
[TF. Following the approach by Molitor et al. [TT, a 
measure for the spatial extent of the localized islands in 
a device can be obtained by finite-bias spectroscopy. In 
devices of different widths, the extracted charging ener- 
gies Ec of the largest diamond in the gap are inversely 
proportional to the width of the nanoribbons and only 
very weakly dependent on the length [l6l[T8]. Again, the 
device investigated in our work falls well into the statis- 
tics of the data published earlier with E^ = 5-10 meV. 

Several microscopic pictures have been introduced to 



explain the formation of a transport gap in graphene 
nanoribbons. Lattice defects at the edges could cause An- 
derson localization |25l [26] which would suppress trans- 
mission around the charge neutrality point. An alter- 
native picture suggests the formation of quantum dots 
along the ribbon due to potential fluctuations [21 . A 
small confinement gap is required in the latter case to pre- 
vent Klein tunneling between the puddles. Experimental 
transport data could so far be interpreted in both mod- 
els. Knowledge about transport mechanisms, which we 
investigate here in thermal activation studies, may help 
to understand where and how localization comes about. 

Fig. [2l[a) displays the back-gate voltage dependent 
conductance at various temperatures for the complete 
transport gap. To obtain this temperature dependence 
of G, the investigated back-gate voltage range was split 
into intervals of about 1 V as indicated in Fig. |2la) by the 
vertical dashed lines. In these sections G was measured 
at stepwise increasing temperatures between 1.25 and 45 
K. In all sections it was verified that the low-temperature 
Coulomb peak spectra were identical before and after the 
thermal cycle. A number of approximately 100 Coulomb 
resonances are visible in the region of suppressed con- 
ductance. For any minimum between two resonances the 
conductance increases for increasing temperatures. Even 
at the highest temperatures the conductance approaches 
but does not exceed e'^/h meaning that the system re- 
mains in the strongly localized regime. 

If we associate with each conductance resonance the 
addition of a single electron to the system, the trans- 
port gap corresponds to a density of states of ^5x10^^ 
m~^eV~^. This value is in good agreement with Ref. 
i35f where the density of states was determined from the 
quantum capacitance of a top-gated large-area device. 

Fig. [2lb) shows a close-up for Coulomb resonances 
with distinctly different behaviors. The amplitude of the 
left peak grows with T and the peak broadens at the 
same time until it is finally swamped away by the rising 
background. The signature of the right peak is a max- 
imum peak value of G at the lowest temperature which 
drops to a local minimum at intermediate temperatures 
and recovers as T is increased further. Such a behavior is 
found only for those ~10 % of the resonances in the inves- 
tigated back-gate window, which are particularly sharp 
at low temperatures. 

For a single quantum dot this observation has been ex- 
plained by the interplay between temperature, the single- 
particle level spacing A^, and the coupling of the en- 
ergy levels to the leads [36 . A strongly coupled ground 
state transition exhibits a 1 /T-dependence for the peak 
height. In contrast, transport through a weakly con- 
nected ground state transition with a strongly coupled 
excitation is enhanced by activation. As kT > A^ both 
levels contribute to transport. 

We now focus on the thermal activation between reso- 
nances. In Fig. [2|c) we display the behavior at three rep- 
resentative back-gate values. In all cases shown here the 
conductance is temperature-independent at low T and 



activation sets in for T > 3 K. The latter is linear in the 
logarithmic plot presented here, which is characteristic 
for activated transport. Thus, the data is fitted to the 
empirical law 



G = Go exp 



E. 



2kBT 



G, 



off, 



(1) 



where Go is the large temperature conductance. Eg, is the 
activation energy and Goff is a constant off-set. Eq. ^ 
is used to fit the T-dependent conductance with those 
three parameters. It reproduces the data very well in 
all conductance valleys between resonances and even on 
some peaks. 
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FIG. 3: Back-gate dependence of fitting parameters: (a) 
Ea, (b) Goff and Go. In (a) and (b) the black solid curve 
shows G at base temperature. Insets: Coulomb diamonds 
reconstructed from Ea for two regimes, (c) Comparison of a 
Coulomb diamond (representing ^c(Vbg)) and E^ determined 
for this VsG-interval. (d) Go{Ea) for the transport gap. Col- 
ored branches indicate Go/E^ pairs that originate from the 
same conductance valleys (arrows in (b)). 



With this model for transport in our system at hand 
we extract Eg,^ Go and Goff as a function of back-gate 
voltage. The analysis was performed only at those gate 
voltages where G{T) spanned more than one order of 
magnitude. The results are shown in Figs, [sj^a) and (b). 
Due to the given criterion for the analysis an evaluation 
at the edge of the transport gap as well as around Fbg 
= -1.25 V was not possible. 



We start with discussing the high-temperature acti- 
vated behavior found in the data. The activation energies 
peak in the middle between neighboring conductance res- 
onances [Fig. [sFa)]. On the other hand, pronounced dips 
in the activation energies arise which coincide with con- 
ductance peaks . In-between a linear dependence on gate 
voltage is observed as is characteristic for Coulomb dia- 
monds. Additionally, the largest Eg, values are 10 to 20 
meV. This energy scale is of the order of typical charging 
energies Ec of this device determined from the finite-bias 
spectroscopy in Fig. [lie). As visualized in Fig. Islc), 
a more careful comparison shows that the activation en- 
ergy resembles the measured Coulomb diamond bound- 
aries remarkably well. Due to thermal cycling in-between 
the diamond and the temperature measurement, some 
shifts are visible in the spectra if the two energy scales 
are plotted on top of each other over a large gate voltage 
range. The finding that the peak values of Eg, in the val- 
leys between conductance resonances are identical to the 
charging energy Ec extracted from Coulomb diamonds is 
a central result of this paper. 

We can reconstruct Coulomb diamonds from the acti- 
vation energy by mirroring Eg,{VBG) at the voltage axis 
and inserting lines along the linear slopes in Eg,. The 
insets of Fig. Isla) displays two qualitatively different re- 
gions in back-gate voltage. In the left graph adjacent di- 
amonds touch each other in one point at zero bias. Their 
size is similar and the flanks have the same slopes. For 
this back-gate voltage range the same observations are 
made for the boundaries of Coulomb blockade diamonds 
measured in finite-bias spectroscopy. Such a behavior is 
characteristic for a single quantum dot where levels are 
filled sequently. In the region under discussion, transport 
is therefore dominated by only one localized island. Since 
the charging happens from the (temperature broadened) 
leads, that are coupled to the island, the corresponding 
maximum Eg, and E^ have to be interpreted as the on- 
site charging energy of this localized site. Its diameter 
corresponds roughly to the ribbon width when estimat- 
ing the size of the puddle from Ec by a comparison to 
data taken on quantum dots. 

The temperature dependence of the conductance res- 
onances between these diamonds exhibits a monotonic 
increase [see Fig. [2la)]. As discussed before, this is ex- 
pected for multilevel transport [36,. 

As a second regime we chose a back-gate voltage in 
Fig. flic) around which the regions of suppressed current 
are connected to each other. The right inset in Fig. [sja) 
shows the corresponding reconstruction of Coulomb dia- 
monds from Eg, where diamonds overlap and the size as 
well as the back-gate dependence of Eg, vary strongly in 
neighboring diamonds. Taking this behavior as an indi- 
cation for the participation of several dots in transport, 
we now have to attribute Eg, and Ec to both on-site and 
inter-site charging energies. Stochastic Coulomb block- 
ade describes such a phenomenon, where transmission 
through a small number of quantum dots is considered. 

Next we proceed with a discussion of the low- 



temperature conductance represented by Goff in eq. ([T]). 
We attribute Goff, which is evident in the curves in Fig. 
[2|c), to cotunnehng processes that determine the conduc- 
tance value before thermal activation sets in. Cotunnel- 
ing leads essentially to Lorentzian tails of conductance 
resonances. The inset of Fig. |2|a) shows that indeed 
we can explain the resonance line shape taking into ac- 
count both thermal and coupling broadening by a con- 
volution of the derivative of the Fermi distribution with 
a Lorentzian. We can do a refined analysis of the low- 
temperature background by fitting the low T data be- 
tween resonances to the expression Giow ^ /3(T^ + Tq ) 
(not shown) ^. Fig. [sj^b) shows that the conductance 
spectrum taken at the lowest temperature is indeed re- 
fiected by the extracted cotunnehng background. The 
finding of cotunnehng transport supports the previous 
statement that only few islands are involved in transport 
since cotunnehng becomes suppressed as the number of 
localized states increases. 

We now discuss the behavior of the prefactor Go in Eq. 
([l]). It extrapolates the conductance for k^T ;» Eg, and 
hence represents the high temperature conductance. The 
order of magnitude of Go is between 0.1 and 1 in units 
of e^ /h. Similar to E^^ it is strongly anti-correlated with 
the conductance at the lowest temperature as illustrated 
in Fig. [sFb). The correlation between Eq^ and Go in 
conductance valleys is visualized in Fig. IsFd). Clearly, 
the Gq{Eq) plot consists of discrete branches with varying 
curvature/slope. Each color-coded branch corresponds to 
a peak of Eg^ in the back-gate spectrum. The ratio of Go 
to £^a decreases as the pair originates from a back-gate 
value closer to the center of the transport gap. 

Transport studies in finite magnetic field have been 
carried out in two devices. As seen in earlier experi- 
ments [17 we find that the size of the Coulomb dia- 
monds shrinks as a 5-field is applied perpendicular to 
the graphene plane. This effect was attributed to time 
reversal symmetry breaking in the regime of strong lo- 
calization which causes an increase of the conductance 
through the ribbon [38^, 'SS^ . To get more evidence for 
the observed positive magnetoconductance we have in- 
vestigated the temperature dependence at 5 = 7 T. As 



for zero field the extracted maximum E^^ is equal to E^ 
of the corresponding Coulomb diamond. 

Comparing the different ribbons under study we ob- 
serve an increase of AFbg and E^ with decreasing rib- 
bon width as discussed in other experiments J9HT9] . The 
latter fact points to the formation of ever smaller islands, 
which block transport and lead to an increase of E^^ as 
the ribbon gets narrower. The magnitude of the intro- 
duced energy £^aC^BG) is tied to E^iV^o) for all measure- 
ments showing that they share the same physical origin. 

Our temperature dependence differs from the one ob- 
served in Refs. ITSl and [TTl where G ex exp(— To/T)^/^ for 
low temperatures. Here, the large number of measured 
points in back-gate allowed us to analyze the temperature 
dependence for discrete V^q values inside the transport 
gap. However, we can fit our data with the same tem- 
perature dependence as in Refs. 15 and 17 if we apply 
the averaging methods used there. 

The picture of transport we present here does not re- 
quire but does not exclude either the contribution of 
phonons inside the system. Activation may take place in 
the leads from which the localized puddles get charged 
via smearing of the Fermi function. It is unclear whether 
phonons in the ribbon get important for transport at 
elevated temperatures. The origin of the correlation be- 
tween £^a and Go remains to be understood but may be 
linked to the role of phonons. 

In summary, we have studied thermally activated 
transport in graphene nanoribbons of different aspect ra- 
tios and compared the determined parameters to trans- 
port measurements at low temperature. We find that 
the transmission is dominated by mainly one of the few 
localized states inside the ribbon at a specific back-gate 
configuration. As a consequence transport in graphene 
nanoribbons should be understood as being mesoscopic 
and single particle-like and treated in such a framework. 
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